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Abstract 

A master equation for the deformed quantum harmonic oscillator interacting with a 
dissipative environment, in particular with a thermal bath, is derived in the microscopic 
model by using perturbation theory, for the case when the interaction is deformed. 
The coefficients of the master equation and of equations of motion for observables 
depend on the deformation function. The steady state solution of the equation for the 
density matrix in the number representation is obtained and shown that it satisfies 
the detailed balance condition. The equilibrium energy of the deformed harmonic 
oscillator, calculated in the approximation of small deformation, does not depend on 
the deformation of interaction operators. 

o ' 



a: 



PACS numbers: 03.65.Bz, 05.30.-d, 05.40.+j, 02.20.Sv 
(a) e-mail address: isar@theory.nipne.ro 



1 Introduction 

> ■ 

For more than a decade a constant interest has been manifested to the study of deforma- 
tions of Lie algebras, so-called quantum algebras or quantum groups, whose rich structure 
produced important results and consequences in statistical mechanics, quantum field theory, 
conformal field theory, quantum and nonlinear optics, nuclear and molecular physics [1, 2]. 
Their use in physics became stronger with the introduction of the g-deformed Heisenberg- 
Weyl algebra (g-deformed quantum harmonic oscillator) by Biedenharn [3] and MacFarlane 
[4] in 1989. There are, at least, two properties which make (/-oscillators interesting objects 
for physics. The first is the fact that they naturally appear as the basic building blocks of 
completely integrable theories. The second concerns the connection between g-deformation 
and nonlinearity. In Refs. [5, 6, 7] it was shown that the g-oscillator leads to nonlinear 
vibrations with a special kind of the dependence of the frequency on the amplitude. The 
g-deformed Bose distribution has been obtained and it was also shown how g-nonlinearity 
produces a correction to the Planck distribution formula [5, 6, 8]. 
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The present paper is the third one in the series of papers devoted to the study of the 
influence of quantum deformation on quantum dissipation. In Ref. [9], using a variant of 
Mancini's model [10], we derived a master equation for the /-deformed oscillator in the 
presence of a dissipative environment for an undeformed interaction between system and 
environment. Then in Ref. [11] we obtained a Lindblad master equation for the ordinary 
harmonic oscillator interacting with an environment through a deformed interaction. The 
equations of motion obtained for different observables had a strong dependence on the 
deformation. In the present paper we set a master equation for the deformed harmonic 
oscillator in the presence of a dissipative environment, for the case of a deformed interaction 
of the system with its environment. This equation is shown to be a deformed version of 
the master equation obtained in the framework of the Lindblad theory for open quantum 
systems [12]. When the deformation becomes zero, we recover the Lindblad master equation 
for the damped harmonic oscillator [13, 14]. We are interested in studying the role of 
nonlinearities which appear in the master equation. This goal is motivated by the fact that 
the g-oscillator can be considered as a physical system with a specific nonlinearity, called 
g-nonlinearity [5, 6]. For a certain choice of the environment coefficients, a master equation 
for the damped deformed oscillator has also been derived by Mancini [10]. 

The paper is organized as follows. In Sec. 2 we remind the basics of the /-deformed 
quantum oscillator, in particular of the g-oscillator. In Sec. 3 we derive a master equation 
for the /-deformed oscillator in the presence of a dissipative environment, for the case of a 
deformed interaction between system and environment. The equations of motion obtained 
for different observables present a strong dependence on the deformation. Then in Sec. 
4 we give the equation for the density matrix in the number representation and find the 
stationary state. In the particular case when the environment is a thermal bath we obtain 
an expression for the equilibrium energy of the oscillator in the approximation of a small 
deformation parameter. A summary and conclusions are given in Sec. 5. 

2 /- and (/-deformed quantum oscillators 

It is known that the ordinary operators {1, a, a', N} form the Lie algebra of the Heisenberg- 
Weyl group and the linear harmonic oscillator can be connected with the generators of the 
Heisenberg-Weyl Lie group. The /-deformed quantum oscillators [15] are defined by the 
algebra generated by the operators {1, A, A', N}, where the Hermitian number operator N 
is not equal to A*>A as in the ordinary case. 
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The /-deformed oscillator operators are given as follows [15]: 

A = af(N) = f(N + l)a, A^ — f(N)at — a}f(N + 1), (1) 
where N = a*a. They satisfy the commutation relations 

[A,N}=A, [A*,N] = -A* (2) 

and 

[A, ^] = (N + l)f(N + 1) - Nf\N). (3) 

The function /, which is a characteristics for the deformation, has a dependence on a defor- 
mation parameter a such that when the deformation disappears, then f(N, a = 0) = 1 and 
the usual algebra is recovered. Transformation (1) of the operators a, a' to A, A* represents a 
nonlinear noncanonical transformation, since it does not preserve the commutation relation, 
i. e. [A,A^]^1. 

The notion of /-oscillators generalizes the notion of g-oscillators (q — 1 + a). Indeed, for 



where q is the deformation parameter (a dimensionless c-number), the operators A and A' 
in Eq. (2) become the g-deformed boson annihilation and creation operators [3, 4]. For the 
g-deformed harmonic oscillator the commutation relation (3) becomes 

[A,A*] = [N + 1]-[N]. (5) 

If q is real positive, then 

= ^Nh3l (6) 
smh(lng) 

and the condition of Hermitian conjugation (Aty = A is satisfied. In the limit q — > 1, 
g-deformed operators tend to the ordinary operators because lim^i [iV] = N. Then Eq. (5) 
goes to the usual boson commutation relation [A, A'] = 1. 

The g-deformed boson operators A and A^ can be expressed in terms of the usual boson 
operators a and (satisfying [a,a^] — 1, N — a) a and [a,N] = a, [at,N] = —at) through 
the relations [16, 17] (see Eqs. (1)): 
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Using the nonlinear map (4) [18, 19], the g-oscillator has been interpreted [5, 6] as a 
nonlinear oscillator with a special type of nonlinearity which classically corresponds to an 
energy dependence of the oscillator frequency. Other nonlinearities can also be introduced 
by making the frequency to depend on other constants of motion, different from energy 
[5, 15], through the deformation function /. Other examples of deformation functions / can 
be found in Refs. [20, 21, 22, 23]. 

The Hamiltonian of the /-deformed harmonic oscillator {uj is the ordinary frequency) is 
a function of N: 

H = ^(AA + A^A) = ((N + 1) f (N + l) + Nf (N)) . (8) 

It is diagonal on the eigenstates \n > in the Fock space and its eigenvalues are 

E n = ^((n+l)f(n + l)+nf(n)). (9) 

In the limit / — > 1 (q — > 1 for g-oscillators), we recover the ordinary expression E n = 
hu(n + 1/2). 

Using the operator Heisenberg equation with Hamiltonian (8) 

dA 

t K— = [A,H] (10) 

we obtain the following solutions to the Heisenberg equations of motion for the /-deformed 
operators A and A^ defined in Eqs. (1) [10, 15]: 

A(t) — exp (—iojQ (N) t) A, A*(t) = A^ exp {iuVt (N) t) , (11) 

where Q(N) is the operator defined as 

Q(N) = X - ((N + 2) f (N + 2) — Nf 2 (TV)) . (12) 

For a g-deformed harmonic oscillator, 

Q(N)= 1 -([N + 2]-[N]) (13) 

and for a small deformation parameter r (r = lng), 

Sl(N) = l + y(iV+l) 2 . (14) 
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3 Quantum Markovian master equation 



In order to discuss the dynamics of the open systems S, we use a microscopic description 
of the composite system. As the subsystem S of interest we take the /-deformed harmonic 
oscillator with Hamiltonian 7i (8) and the environment R (reservoir, bath) with the Hamil- 
tonian Hr. The coupled system S+R with the total Hamiltonian = 7i + Hr + V (V 
is the /-deformed interaction Hamiltonian) is described by a density operator which 
evolves in time according to the von Neumann-Liouville equation 

!<gU-i[tf T ,xC)l. (15) 

When the Hamiltonian evolution of the total system is projected onto the space of the har- 
monic oscillator, the reduced density operator of the subsystem is given by pit) = TrRx(t). 
The derivation of the reduced density operator in which the operators of the environment 
system have been eliminated up to second order of the perturbation theory can be taken 
from literature [24, 25, 26, 27, 28]. Following [27] and the procedure used in the previous 
paper [9] , we point the main steps in obtaining the master equation which describes the time 
evolution of the damped deformed harmonic oscillator. We assume that the initial state i?o 
of the environment at t — does not depend on the state p(0) of the subsystem. Then 
x(0) = p(0)R . At later times correlations between S and R arise due to the coupling of 
the system and environment through V. However, we assume that this coupling is weak and 
then at any instant of time, not only for the initial time, the reservoir and the oscillator are 
approximately decoupled. Furthermore, R is a large system whose state should be virtually 
unaffected by its coupling to S. We then write x(t) = p(t)R . In the previous paper [9], 
we considered the interaction potential of the linear form in the coordinate and momentum 
operators in the Hilbert space of the subsystem S. In analogy to that model, we assume in 
the present paper a /-deformed interaction potential V of the form 

V= £ SiFi, (16) 

i=l,2 

where S± — Q and S2 = P are the /-deformed coordinate and momentum operators in the 
Hilbert space of the subsystem S (m is the oscillator mass): 

with A', A defined in Eqs. (1) and Tj are Hermitian operators in the Hilbert space of the 
environment. Then we obtain the following master equation for the density operator of the 
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open quantum system in the Born-Markov approximation: 
dp(t) 



dt 



-Un, P (t)] + i E tdt'ic^t'^pmi-t')} + cMs j (-if)p(t),s i Mi8) 



where the coefficients Cij(t') = Tr ji{RoTi(t')Tj} are correlation functions of the environment 
operators. It is assumed that the correlation functions decay very rapidly on the time scale 
on which p(t) varies. Ideally, we might take Cij(t') ~ S(t'). The Markov approximation relies 
on the existence of two widely separated time scales: a slow time scale for the dynamics 
of the system S and a fast time scale characterizing the decay of environment correlation 
functions [27]. 

In order to get the time dependence of the operators of coordinate S\(t) = Q(t) and 
momentum S 2 (t) = P(t) used in Eq. (18), we express them through the relations 



Q(t) = ^^(A\t) + A(t)), P(t) = ij^(A\t)-A(t)) (19) 

and then insert Eqs. (11) for A(t) and A^(t). Then the master equation (18) takes the 
following form: 

+i f dt'(C* n (t')[Q, pit) (qE^ + E+Q- — {PE_ - E+P))] 
2n Jo V muj J 

+C n (t')[(QE_ + E+Q- — (PE. - E+P)) p(t), Q] 
+iC* 2 (t')[P, p(t)(mw(QE_ - E + Q) - i{PE_ + E+P))} 
+iC 22 {t')[{mu (QE- - E + Q) - i {PE„ + E + P))p(t),P] 
+iCl 2 {t')[Q,p(t){mw (Q£_ - E + Q) - i (PE_ + E + P))} 

+C 2 i(t')[(QE- + E+Q--^- (PE_ - E+P)) p(t), P] 

+C* 21 (t')[P,p(t) (qE_ + E+Q- — (PE_ - E+P))} 
+iC 12 (t')[(mu {QE. - E+Q) - i (P£_ + E+P))p(t),Q}), (20) 

where we have introduced the following notations: 

E+ = exp(iujn(N)t'), E- — exp (—iuQ (N) t') . (21) 

If the environment is sufficiently large we may assume that the time correlation functions 
decay fast enough to zero for times t' longer than the relaxation time ts of the environment: 
t' >> ts- Therefore, if we are interested in the dynamics of the subsystem over times which are 
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longer than the environment relaxation time, t ^> we may use the Markov approximation 
and replace the upper limit of integration t by oo. Physically, this amounts to assuming 
that the memory functions Cij(t')E±(t') decay over a time which is much shorter than the 
characteristic evolution time of the system of interest. We evaluate the integrals in Eq. 
(20) like in the preceding paper [9]. After certain assumptions [24, 28], one can define the 
complex decay rates, which govern the rate of relaxation of the system density operator, as 
follows: 



In fact, D PP (Q), D qq (Q), D pq (Vt) and X(Q) play the role of deformed diffusion and, respec- 
tively, dissipation coefficients. The existence of these coefficients reflects the fact that, due 
to the interaction, the energy of the system is dissipated into the environment, but noise 
arises also (in particular, thermal noise), since the environment also distributes some of 
its energy back to the system. Then the master equation (20) for the damped deformed 
harmonic oscillator takes the form: 



Here the curly parentheses mean anticommutators. We notice that the deformation is 
present in both the commutator containing the oscillator Hamiltonian Ti, as well as in 
the dissipative part of the master equation, which describes the influence of the environ- 
ment on the deformed oscillator. This master equation preserves the Hermiticity property 
of the density operator and the normalization (unit trace) at all times, if at the initial time 
it has these properties. In the limit / — > 1 (fl — > 1), the deformation disappears and Eq. 
(25) becomes the Markovian master equation for the damped harmonic oscillator, obtained 
in the Lindblad theory for open quantum systems, based on completely positive dynamical 




(22) 




(23) 




(24) 
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semigroups [12, 13, 14]. The fact that we introduced a /-deformed interaction between the 
open sytem and environment is reflected in the presence of deformed operators P and Q 
in the dissipative part of Eqs. (20) and (25), along with the deformed diffusion operator 
coefficients, compared to the previous model [9, 10], where the dissipative part of the master 
equation contains the deformed diffusion coefficients and the usual operators p and q. 

Expressing the coordinate and momentum operators back in terms of the creation and 
annihilation operators, introducing the notations 

D+la) . i_ ^ Dm (n) + , D _ (a) s _L (n) _ (26) 

and assuming, like in Ref. [9], that A(f2) = A = const, the master equation (25) for the 
damped deformed harmonic oscillator takes the form 

f = ~[n, P ] + ([[D + (Q)af(N),p}J(N)a^} 

-[[f(N)a* (D_ (fl) + % -D pq (fl)) , p], f{N)a)] - ^[f{N)a\ {af(N),p}} + H.c), (27) 

with H given by Eq. (8). 

In the particular case of a thermal equilibrium of the bath at temperature T (k is the 
Boltzmann constant), we take the diffusion coefficients of the form (in concordance with 
results of Refs. [9, 10]) 

mw D w (n) = ^^ = |Aooth^, D pq (Q) = (28) 
and the master equation (27) takes the form 

f = ~[n, P ] + ± {[[coth^af(N),p}J(N)a^} - \f(N)a\ {af(N), p}} + H.c^j . (29) 

In the limit Q — > 1, the deformed diffusion coefficients (28) take the known form obtained 
for the damped harmonic oscillator, if the asymptotic state is a Gibbs state [13, 14]: 

Timuj huj h Tiu) 
Dvv — Acoth— — , D aa = Acoth— — -, D vo = 0. (30) 

If the bath temperature is T = 0, the master equation (29) simplifies: 

f = -y l(N + l)f(N + l) + Nf(N),p] 
-A (Nf 2 (N)p + pNf(N) - 2f(N + l)apa) f(N + 1)) . (31) 

The meaning of the master equation (29) becomes clear when we transform it into 
equations satisfied by the expectation values of observables involved in the master equation, 
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< O >= Tr[p(t)0], where O is the operator corresponding to such an observable. We give 
an example, multiplying both sides of Eq. (29) by the number operator iV and taking the 
trace. Then we obtain the following equation of motion for the expectation value of N: 

j t <N>= A(< (coth - 1)(AT + l)f 2 (N + 1) > 

- < (eot. /'^- 1 W / 2 W >). (32) 

This equation leads to a time dependence of the averaged number of quanta on dissipation, 
temperature and deformation, compared to the case of an oscillator without dissipation, 
where the expectation number of quanta is conserved. In the case of a g-deformation, Eq. 
(32) takes the form 

| < N >= A (< (coth -!)[*+!]>-< (coth MM^il + 1)[N] >) . (33) 

If, in addition, the temperature of the thermal bath is T = 0, Eq. (33) becomes 

£ < N >= -2A < [N] > . (34) 

In order to obtain an approximate solution of this equation, we work in the limit of a small 
deformation parameter r = lng (q real). Then we can take [1, 2] 

[N] = N + — (N 3 -N) (35) 

and, making also the assumption < iV 3 >m< N > 3 , Eq. (34) reduces to the following 
differential equation: 

A < N >= -2A ( < N > +^(< N > 3 - < N >) ) . (36) 
at \ 6 ) 

By integrating this equation we obtain 

< N(t) > < N(0) > e- 2 ^ 1 -^ 



V /l -T + T< Ar W > 2 V /l -T + T< Ar (°) > 2 



(37) 



from where we can obtain the following expression for the expectation value of the number 
operator in the approximation of a small deformation parameter: 



< N(0) > Jl - 2± e -2A(l-V)* 

< N(t) >= 



1 - T + T < iV (°) >2 ( X ~ e-^-V)*) 
t< N(0) > e- 2 ^ 1 -^' (l~^< N(0) > 2 (1 - e- 4 ^ 1 - 2 ^')) • (38) 
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In the limit r — > 0, we obtain 

< N(t) >=< N(0) > e- 2Xt , (39) 

which is the expression of the expectation value of the number operator of the damped 
harmonic oscillator obtained in the Lindblad theory for open quantum systems. 

We consider another example, taking the simplest case of a thermal bath at T = 0, 
when both diffusion and dissipation coefficients do not depend on the deformation, D + = 
A/2 = const. Even in this situation, the equations of motion for the expectation values are 
yet complicated, because they do not form a closed system. Multiplying both sides of Eq. 
(31) by the operator a and, respectively, fl(N)a and taking throughout the trace, we get 
the following equations for the expectation values of these operators: 

-f- < a >= —iu < Q(N)a > 
at 

-A < ((N + l)f 2 (N + 1) + Nf 2 (N) - 2Nf(N)f(N + 1)) a >, (40) 

4 < £l(N)a >= -iu < n 2 (N)a > 
at 

-A < (n(N) [(N + 1) f (N + 1) + Nf (TV)) - ^{N - l)Nf(N)f(N + 1)) a > . (41) 



These examples show that the equations of motion contain strong nonlinearities introduced 
by the deformation function / and they do not form a closed system of equations. 

4 Equation for the density matrix 

Let us rewrite the master equation (27) for the density matrix by means of the number 
representation. Namely, we take the matrix elements of each term between different number 



states denoted by \n >, and using N\n >= n\n >, a + \n >= \/n + l\n + 1 > and a\n >- 
y/n\n — 1 >, we get 

= - yLWV) + (rn + l)/ 2 (m + 1) - nf 2 {n) - (n + l)f{n + l)\p mn 
-[(m + l)f 2 (m + 1) (D + (Q(m)) - £j + mf(m) (D + (n(m - 1)) + ^ 

+(n + l)f 2 (n + 1) (i? + (n(n)) - ^ + nf 2 (n) ( D + (Q(n - 1)) + £ )}p mn 



+yj(m + l)(n + l)/(m + l)/(n + l)[D+(ft(m)) + D+(fi(n)) + A]p m+1 , n+1 
+v^/(m)/(n)[£)+(fi(m - 1)) + D+(fi(n - 1)) - A]p m _ lin _! 
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■yj(m + l)nf(m + l)/ 2 (n)[D_(fi(m)) + - 1)) 



1 

-- (D pq (tt(m)) + D pq (Q(n - l)))]p m+ i, n _i 



-yjm(n + i)f(m)f(n + l)[LL(fi(m - 1)) + £>_(ft(n)) 
+ - (£> M (ft(m - 1)) + £> M (ft(n)))]p 

m— l,n+l 

+ v /(m + l)(m + 2)/(m + l)/(m + 2)[ J D_(fi(m + 1)) - ^D pq (Q(m + l))]p m+2 , n 
+ v /(n + l)(n + 2)/(n + l)/(ra + 2)[£L(ft(n + 1)) + ^ OT (fi(n + l))]p m ,„ +2 
+ v /m(m-l)/(m - l)/(m)[£L(ft(m - 2)) + l -D pq (Q(m - 2))]p m _ 2 , n 

+^Jn(n - l)f(n - l)/(ra)[£L(ft(n - 2)) - ^ M (H(n - 2))]p m , n _ 2 . (42) 

Here we used the abbreviated notation p mn =< m\p{t)\n > . This equation, very complicated 
in form and in indices involved, gives an infinite hierarchy of coupled equations for the matrix 
elements. When 

£L(n(n)) = 0, D pq (Q(n)) = 0, (43) 

the diagonal elements are coupled only amongst themselves and not coupled to the off- 
diagonal elements. In this case the diagonal elements (populations) satisfy a simpler set of 
master equations: 

¥M = - (( n + l)f{n + 1) (2£> + (n(n)) - A) + nf{n) (2D + (Q(n - 1)) + A)) P(n) 
+ (n + l)f{n + 1) (2D + (n(n)) + A) P(n + 1) + nf\n) (2D + (Q{n - 1)) - A) P(n - 1), (44) 

where we have set P(n) = p„ n . For a g-deformation, Eq. (44) takes the form 

= _ ([„ + 1] (2£> + (n(n)) - A) + [n] (2L> + (fi(n - 1)) + A)) P(n) 
+ [n + 1] (2D+(fi(n)) + A) P(n + 1) + [n] (2D + (Q(n - 1)) - A) P{n - 1), (45) 



where, according to Eq. (13), 



n(n) = i([n + 2]-[n]). (46) 



We define the transition rates 

*+(n) = (n + l)/ 2 (n + 1) (2£>+(ft(n)) - A) , t_(n) = n/ 2 (n) (2£>+(ft(n - 1)) + A) , (47) 
which for a g-deformation look like 

t + (n) — [n + l] (2D + (Q(n)) — A) , f_(n) = [ra] (2Z>+(ft(n - 1)) + A) . (48) 
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With these notations Eq. (44) becomes: 

d ^^- = t + {n - l)P(n -l) + t-(n+ l)P(n + 1) - (t+(n) + tJn)) Pin). (49) 
at 

The steady state solution of Eq. (49) is found to be 

We remark that in the steady state the detailed balance condition holds: 

t-(n)P(n)=t+(n-l)P(n-l). (51) 

In the particular case of a thermal state, when the diffusion coefficients have the form (28), 
the stationary solution of Eq. (49) takes the following form: 

Pfsin) = Zf exp{-^ ((„ + l)f(n + 1) + nf(n)) }, (52) 

where 

ZJ^PMexpM^l (53) 
and Zf is the partition function: 

oo * 

Zf = T, ex P{"^ + l )^ n + !) + ^ 2 H)>- ( 54 ) 

By using the eigenvalues (9), the distribution (52) can be written 

^(n) = Z / - 1 exp(-^). (55) 

Expressions (52) and (55) represent the Boltzmann distribution for the deformed har- 
monic oscillator. In the limit / — > 1 the probability -P*^(n) becomes the usual Boltzmann 
distribution for the ordinary harmonic oscillator with the well-known partition function 
Z = 1/2 smh^. For the g-oscillator described by Hamiltonian H in Eq. (8) and weakly 
coupled to a reservoir kept at the temperature T, the (/-deformed partition function can be 
obtained as a particular case of the partition function Zf (54), by taking the deformation 
function (4): 

Z ^gexp{-|^ sinh(r( '' +1) » +sinh( "' ) }. (56) 
n _Q Zrci sinn T 

The results (50) - (56) coincide with those obtained in the previous paper [9], where we 
considered a linear interaction potential V and the harmonic oscillator operators contained 
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in this potential are kept undeformed. Therefore, for the equilibrium energy we obtain in 
the limit of small deformation parameter r the same expression like that obtained in Ref. 
[9]: 

E(t - oo) = — (coth — + r 2 c), (57) 

where 

c= (^iy(iF3i-^ {eP _ l? J. = ^- ( 58 ) 

We note that in the approximation of small deformation, the energy of the damped deformed 
oscillator depends on the energy hu/2 of oscillator ground state and on the temperature T. 
Evidently, when there is no deformation (r — > 0), one recovers the energy of the ordinary 
harmonic oscillator in a thermal bath [13, 14]. In the limit T — > 0, one has c — > 0, E(oo) = 
Tiuj/2 and the deformation does not play any role. 



5 Summary and conclusions 

Our purpose was to study the dynamics of the deformed quantum harmonic oscillator in 
a deformed interaction with a dissipative environment, in particular with a thermal bath. 
We derived in the Born-Markov approximation a master equation for the reduced density 
operator of the damped /-deformed oscillator. The one-dimensional /- or g-oscillator is 
a nonlinear quantum oscillator with a specific type of nonlinearity and, consequently, the 
diffusion and dissipation coefficients which model the influence of the environment on the 
deformed oscillator strongly depend on the introduced nonlinearities. Compared to the 
previous paper [9], in the present work the harmonic oscillator operators contained in the 
interaction terms are deformed. This fact is reflected in the increased degree of the non- 
linearity of the equations of motion for the expectation values of observables. In the limit 
of zero deformation, the master equation takes the form of the master equation for the 
damped oscillator obtained in the framework of the Lindblad theory of open systems based 
on quantum dynamical semigroups. We have also derived the equation for the density ma- 
trix in the number representation. In the case of a thermal bath we obtained the stationary 
solution, which is the Boltzmann distribution for the deformed harmonic oscillator. In the 
steady state the detailed balance condition holds true. In the approximation of a small de- 
formation, we obtained the expression of the equilibrium energy of the deformed harmonic 
oscillator without deformation of interaction. This energy depends on the oscillator ground 
state energy and temperature. 
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The master equation for the damped deformed harmonic oscillator is an operator equa- 
tion. It could be useful to study its consequences for the density operator by transforming 
this equation into more familiar forms, such as the partial differential equations of Fokker- 
Planck type for the Glauber, antinormal ordering and Wigner quasiprobability distributions 
or for analogous deformed quasiprobabilities [7] associated with the density operator. It 
could also be interesting to study the properties of the entropy of the damped deformed 
harmonic oscillator and to find the states which minimize the rate of entropy production for 
this system. In the case of the undeformed damped oscillator such states are represented 
by correlated coherent states [29] . We suppose that for the damped deformed oscillator the 
corresponding states are the deformed (nonlinear) coherent states [15], which play an impor- 
tant role in the description of the phenomenon of environment induced decoherence. The 
dissipative dynamics of deformed coherent states superposition and the related coherence 
properties have recently been studied by Mancini and Man'ko [30]. 
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